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■ , Abstract 

CIh' The baryon number dissipation rate due to sphaleron transitions at high 

r~| , temperatures in the minimal standard model is evaluated. We find that this 

rate can be considerably suppressed by one loop contributions of bosonic and 
^ . fermionic fluctuations which are particularly important for a small mass of the 

j^ ■ Higgs boson and a large top quark mass. Fixing the latter to its recently stated 

value of 174 GeV the complete erasure of the baryon asymmetry is prevented 

within the framework of the minimal standard model if the Higgs mass is less 

than about 66 GeV. 
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1. The potential energy of the SU{2) Yang-Mills gauge field and the Higgs field 
in the minimal standard model is periodic in the Chern-Simons number Nqs 
[ffl, 0|. The vacua (with integer Ncs) are separated by the sphaleron barrier 
[^, Q whose height is of the order of mw/a, where niw is the VF-boson mass 
and a = g^ /{At:) is the SU{2) gauge coupling constant. 

Transitions from one vacuum to a topologically distinct one over this barrier 
cause a change in the baryon and lepton number by one unit per fermion 
family due to the axial anomaly |^. Hence, this transition is a baryon and 
lepton number violating process. Although it is very strongly suppressed under 
ordinary conditions |5[, its rate can become large at high densities [^, 0], high 
temperatures ||8|, ^, IC] or maybe, at high particle energies [|l^]. In particular. 



the transitions should have occurred in the early universe immediately after 
the electroweak phase transition. Whatever was the mechanism leading to the 
baryon asymmetry of the universe at earlier times, these transitions might have 
washed out any initial excess of baryons over antibaryons if their rate was large 
enough |^, ^, |l2[ . In order to understand the presently observed baryon excess 
we therefore need an exact determination of the transition rate. 

The dominant contribution to this rate is given by the classical Boltzmann 
factor e~ '='''==' where T is the temperature, and -Eciass is the classical energy 
of the sphaleron. Quantum corrections arise from bosonic and fermionic fluctu- 
ations about the sphaleron; there are also prefactors due to negative and zero 
bosonic fluctuation modes |]l^]. In |1^ the rate was calculated considering the 
classical contribution and the factors due to zero and negative modes only; in 
this case the rate is so large that any initial baryon excess could easily have 
been washed out. Hence it is necessary to take the loop corrections into ac- 
count if one wants to have a chance to preserve the asymmetry. In 1^, ^ boson 
fluctuations were considered through an effective Higgs potential resulting in 
an upper limit of 45 to 55 GeV for mn- 

A direct evaluation of the boson determinant over non-zero modes was made 
in [13 1 where an approximation method |14] was employed; exact calculations 



were performed in |15, 16 1. All those calculations were done in the limit of 
high temperature in which the four dimensional boson fluctuation matrix can 
be replaced by the three dimensional one while fermions decouple completely. 
Although parametrically this limit is reasonable, numerically it might not nec- 
essarily be justified. In fact we find that numerical results are significantly 
influenced by terms which vanish in the high T limit, especially by the contri- 



bution of the fermion fluctuations. Moreover, a third independent calculation 



of the boson determinant is necessary since the results of |15l and [16] deviate 
from each other. 

We have consistently evaluated both the boson and the fermion one loop 
contributions by the determination of the complete (discretized) spectrum of 
the fluctuation matrices. This technique allows us not only to check the existing 
results in the high T limit but also to perform the generalization to arbitrary 
temperatures. A detailed description of this technique for the fermions can be 
found in M, an extended paper about the boson fluctuations is in preparation 



2. We consider the minimal version of the standard electroweak theory with 
one Higgs doublet which is Yukawa coupled to left handed fermion doublets and 
to right handed singlets; in the following we write only one doublet and one 
pair of singlets for brevity. We neglect the Weinberg angle, i.e. we work with 
a pure SU{2) gauge theory. This idealization does not seem to be significant 
[p^]. The Lagrangian is thus 

+ i^dl^D^iljL + XRij^df.XR - i^LMXR - XrM^l (1) 

with F" = df^A"^ — d,^A°j^ + e"'^'^ A^ Af, and the covariant derivative being defined 
as D^ = d^—^A'^r"'. Af is a 2x2 matrix composed of the Higgs field components 
$ = (^o) and the Yukawa couplings /i„, /i^: 

M=( '"•*°' "'^M. (2) 

11)1 means the SU{2) fermion doublet (^^), and with Xr we denote the pair of 
the singlets ur, (Ir. To preserve spherical symmetry we take equal masses for 
the two kinds of fermions in each doublet, i.e. hu = hd = h. This approximation 
is justified for the light doublets, but not so good for the bottom and top quark 
doublets. Here the mass difference could be treated as a perturbation. However, 
preliminary estimates show that the final conclusions of the present paper do 
not change noticably. 

The transition rate per volume V of the system going from one vacuum 
into a topologically distinct one is given by the semi-classical Langer-Affleck 



formula [^, 2C, ^ which, apphed to the model ||T^, reads 

r 



_ _uj- A/o - 

V 111 V 



exp(-^class/7') , 



(3) 



with uj- being the frequency of the negative bosonic mode and Nq- the volume 
factor due to the negative and the zero bosonic modes. -Eciass is the classical 
energy of the sphaleron configuration, 



ii^class = jd\ [-L(F^.)' + (A^)^(A^) + y ($^$ 
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The determinants Kbosi '^fcrm correspond to bosonic and fermionic fluctuations; 
they are divergent and have to be renormalized which we perform using a 
proper time regularization scheme. The divergent parts are combined with the 
classical energy, which yields the physical parameters normalized at the scale 
of mw- The major part of their temperature dependent contribution can also 
be absorbed by the classical energy, leading to temperature dependent masses 
(see e.g. [pi): 



mw{T) mniT) mpiT) Ecii,ss{T) 



QiT) 



'1 



2^2 

if' 



mw ruH mp -Edass 

where the critical temperature of the phase transition is given by [^, |^ 
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In a realistic model all fermion masses except the top mass are negligible, so 

2 
doubl. ''^F 



that X^doubi. ''^F ~ i"^t • After performing this renormalization we find 
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with A'^rot, A'trans being the Jacobians of the zero modes. The determinants can 
be written in the following form: 
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Here e^ are the eigenvalues of the fermionic fluctuation operator S'^Scg/dipSip 
[0], u!^ are the eigenvalues of the bosonic fluctuation operator, using the R^=i 
background gauge for the fluctuations |13, 16|, and (w^^)^ are the correspond- 
ing eigenvalues of the Faddeev-Popov operator. For numerical purposes the 
spectra are discretized and made finite using suitable box parameters which 
have been taken large enough to ensure stability. The double prime Yl", J2" 
means that negative and zero mode frequencies are removed. Knowing the 
eigenvalues e^, cOn and co^^ the above determinants eqs. (^, |9|) can be evaluated 
for any temperature T. The first terms in the exponents correspond to the 
zero temperature fiuctuation energy while the logarithms are only present for 
non zero temperature. In previous works |13, ^, 16 1 only the high T limit was 




considered, in which the determinants reduce to 

'^ferni 

'^bos = ( -^ I ( ^^^ I ( ;^r^Sk7T I • (10) 

We assume the hedgehog ansatz 

.a. X 1-A{r) ^^ R(r) C{r) 

Ai (r) = EaijUj h [Oai - naUi) h UaUi , 

cl>(r) = -|[i/(r)+iG(r)n.0]O, (11) 

for the classical fields and solve the corresponding equations of motion for 
the profile functions A{r), . . . ,G{r). Then the fluctuation operators become 
block diagonal with respect to some "grand spin" K and can be numerically 
diagonalized for each K separately. More details can be found in |Q, |l^]. Hence 
it is possible to evaluate the rate 7(r) numerically for given values of mn /mw 
and mt/m\Y as a function of the temperature. 

Before we present numerical results we would like to get a feeling for the 
dependence of the determinants on the mass parameters. Though the exact 
numerics is quite complicated, it is easy to find the behaviour of the determi- 
nants for large top quark mass and for small Higgs mass. For large top quark 



mass the square loop diagram in the external Higgs field is the dominant con- 
tribution, and the fermion sea energy is proportional to Nc{h(j))^ln{h(j)/mw), 
where h is the Yukawa coupling and (j) is the Higgs field of the sphaleron. To 
get the energy one has to integrate this expression over the range of space 
where the Higgs field differs from its vacuum expectation value, i.e. over the 
spread of the sphaleron. For large spatial distance r the asymptotics of the 
fields is dominated by the term e~^"^ so that the sphaleron size is roughly 
~ m^ if TTiH < mw- Hence all matrix elements of the fiuctuation matrices 
behave as mjj . This scaling law holds for all contributions to the determinant; 
i.e. for small Higgs masses the fluctuations show a strong increase. It is this 
increase which is finally responsible for the suppression of the baryon number 
dissipation and provides an upper limit for niH to ensure the preservation of 
the asymmetry. Thus for mn < rnw < i^t we can estimate 

ln(K;fcrm) ~ Nc z -. r, , (12) 

{mH/mwr 

These are the quantum corrections to the classical sphaleron energy Sciass 0- 
Though parametrically they are a times smaller than i^ciass i numerically the 
fermion sea contribution to the sphaleron energy appears to be large, especially 
for large top masses mt and relatively small Higgs masses mn- It is mainly the 
fermionic factor Kferm j which was put to unity in the previous work P, ^, 15, 



p^ ] , that leads to a significant additional suppression of the baryon dissipation 
rate, see below. 

3. Sphaleron transitions can increase and decrease the baryon number. If the 
baryon number B were zero, the transitions in both directions would happen 
equally often and cancel each other. In the case B ^ one has to introduce a 
chemical potential; it favours transitions which erase the baryon asymmetry in 
accordance with the le Chatelier principle. This has been done in |9|, ^], and 
for fermions with small masses one gets: 

]_dB^_13j(Tl 

B dt 2 r3 ■ ^ ^ 

Since the top mass is actually not small the prefactor ^ should be replaced by 
a slightly bigger number. This effect, however, is negligible compared to the 
other factors so we do not consider it further. 
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In this letter we assume, in accordance with the standard model, that B — L 
is conserved; we do not consider a possibility that there might be a primodial 
excess of say, antileptons owing to unknown forces violating the B — L number. 
In that case the sphaleron transitions would, on the contrary, lead to the yield 
of baryons. 



Standard cosmology gives a relation between time and temperature |23|: 

t = CT-^, (15) 

with the constant C ~ 5 • 10^^ m]y depending on the Planck mass and the 
number of degrees of freedom of the thermalized particles. Hence we obtain 

J_^ = 13C31P, (16) 

B dT r6 ' ^ ^ 

which can be integrated to 



B{T) = B{Tc) e^p\ -13 cJ^^^dT 



where q{T) was defined in eq. (|5|). Thus, we can evaluate the ratio B{0)/B{Tc) 
of the present baryon number, B{0), to the one immediately after the phase 
transition, B{Tc). From this ratio we shall conclude on an upper bound of the 
Higgs mass. 

4. For our numerical calculations we fixed the coupling constant to its physical 
value g = 0.67, mw = 83 GeV, and the top quark mass to mt = 2.1 mw , 
according to its recently stated value of 174 GeV. The only free parameter left 
is the Higgs mass mn- 

First we recalculated In(Kbos) of eq. ( |lO| ) in the high T limit which was 
previously done in |15, Il6| with deviating results. We find that our values 



are between those of the two preceeding computations; they agree with the 
ones of ref. [|^ up to about 10%. Apart from numerical uncertainties this 
difference could be probably explained by the fact that in ref. |l^ the complete 
tadpole expansion except the term linear in temperature was removed due 
to the loop renormalization, while we have performed the renormalization at 
zero temperatures strictly, as it is usually done. Our scheme corresponds to a 
subtraction of only the first term in the tadpole expansion so that the difference 



lies in higher order terms which vanish in the high T hniit but can influence 



the numerical result. There is a larger deviation from the results of [15|, only 
a qualitative agreement is found. 

A powerful check of our numerical performance was that determinants were 
computed in several completely different gauges for the sphaleron field. The 
fluctuation operators change drastically when one goes from on gauge to an- 
other, however we have checked that the gauge invariant eigenvalues change 
only in the range of 0.5%, giving an estimate of the numerical accuracy of the 
spectrum. 

Our aim is, however, a generalization to arbitrary temperatures. Figure 1 
shows the logarithms of the different contributions to the transition rate 7 of 
eq. @: the classical part e~^<=i'»==(-^)/-^, the prefactor uj^N'o-/{2-nV), and the 
bosonic and fermionic fluctuations Kbos and Kferm- It is convenient to take the 
parameter q{T) (see eq. (^) as independent variable rather than the tempera- 
ture itself. Obviously, the loop corrections, especially the fermionic ones yield 
a strong suppression of the total transition rate. If the temperature approaches 
the critical value of the electroweak phase transition, all masses disappear, 
and the transition rate goes rapidely to zero due to the vanishing prefactor. 
It should be mentioned that the one-loop approximation we are dealing with 
breaks down in the near vicinity of the phase transition. However, the integral 
in (17) is strongly dominated by a small region around the maximum of the 
transition rate 7(9), which is separated from the phase transition itself (see 
Figure 1). Figure 1 was calculated for niH = Q.Smw, for other values of mn 
the qualitative behaviour of the curves is basically the same, only the numbers 
change somewhat. 

According to eq. ( |T7| ) we calculated the quotient Bq/Bt^ of the present 
baryon number Bq and the initial one Bt^ . Figure 2 shows the ratio for different 
values of mn ■ If we assume that the initial baryon number Bt^ was not larger 
than about W^ Bq [^], we find an upper bound of the Higgs mass of about 
0.79 mw ~ 66GeV. Figure 2 exhibits a very strong dependence on rriH'- even 
if we assume an apparently unrealistic initial baryon excess Bt^ ~ 1, we get 
from the present-day asymmetry Bq ~ 10^^*^ that the Higgs mass does not 
exceed the value of niH ~ 0.81 mw ~ 67 GeV. On the other hand, if we assume 
Bt^ = Bq, the upper limit for mn lies at mn ~ 0.7mw ps 58 GeV. 

In order to emphasize the significance of the fermionic fluctuations, which 
have not been taken into account previously, we calculated the transition rate 



and the resulting decrease of the baryon number also without the ferniionic 
contribution. Figure 2 shows that in this case the upper bound for the Higgs 
mass, assuming Bt^ = lO^i^O; would be as low as about 51 GeV which is close 
to the approximate estimate of ref. y, 12]. 



To summarize: For any temperature T < Tc both bosonic and ferniionic 
fluctuations suppress the sphaleron transition rate considerably, depending on 
the value of the Higgs mass. Setting the top quark mass to its recently claimed 
value of 174 GeV, the condition that the baryon asymmetry is not washed 
out by sphaleron transitions leads to an upper limit of 66 GeV for the Higgs 
boson mass. This numerical result is based on the Langer-Affleck formula and 
applies to one loop calculations within the minimal standard model with one 
Higgs doublet. Our upper limit is compatible with the experimental lower limit 
of about 60 GeV |24, |2^ provided we include the fermionic corrections. If no 



Higgs particle is found below our upper limit, this would be a clear sign that 
the minimal standard model should be extended, maybe to a model with two 
Higgs doublets. 
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Figure Captions 

Figure 1: Classical and loop contributions to the sphaleron transition rate per 
volume, q = \/l — {T /Tc)"^ , Tc = 0.957 mw, rriH = 0.8 mw ■ The interval 
which actually contributes to the integral ([l7|) is marked. 

Figure 2: Total decrease of the baryon number via the sphaleron transition in 
dependence of the Higgs mass, with and without fermionic fluctuations taken 
into account. Bq/Bt^ is the ratio of the present baryon number to the one of 
the electroweak phase transition. 
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Contributions to transition rate 
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